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1 A description is given of an apparatus used for determining the thermal conductivities 
of loose-fill granular materials in the temperature range 100 to 1,100 °C. The method used 
in making these determinations was that of radial heat flow through a hollow cylinder of 
specimen material contained between a central ceramic core and an outer ceramic shell. 
The heat flow through the specimen was determined by measurement of the power input to 
a heater in the central ceramic core. The radial temperature drop through the specimen 
was inferred from temperatures measured in the core and in the shell, thus avoiding entirely 
the problems of measuring a tcunperature difference within the specimen material. Experi- 
mental results with an estimated uncertainty of 8 percent or less are presented for diato- 
maceous earth having a density of 0.15 g/cm^, and for aluminum oxide powders having 
densities of 0.40 and 0.44 g/cm^. A discussion of the test method is given, with attention 
to possible sources of error, f 



1. Introduction 

Various loose-fill thermal insulations are used, in 
conjunction with appropriate t>uar(lin<>-, to reduce 
extraneous heat flows in several thernuil conductivity 
uieasurino- devices under development in the Heat 
Transfer Section at the National Bureau of Stand- 
ards. Shice such heat flows cannot h(* entirely 
eliminated, it is necessary to know the ai)proxiniate 
thermal conductivity of these insulations in ov(\vv to 
effect appropriate corrections. In designinii' an 
apparatus Tor determining the thermal conductivity 
of these insulations, it was decided to aim for an 
accuracy of 5 percent or better, because of the 
general interest in thermal conductivity of loose-fill 
materials, both from an engineering and from a 
theoretical viewpoint. 

The purpose of this paper is to describe an appa- 
ratus and method used for determining the thermal 
conductivity of loose-fill materials. Experimental 
results are presented for diatomaceous earth and 
alumirmm oxide powder in the temperature range 
from 100 to 1,100 °C. Factors affecting both the 
precision and accuracy of the measurements are 
discussed. 

2. Method 

Radial heat flow in a hollow cylinder of loose-fill 
insulating material was selected as the method for the 
measurements. A cross section of the essential 
elements of the apparatus is shown in figure 1 . The 
specimen was contained within the annular space 
between the outer radius, a, of the ceramic core and 
the inner radius, b, of the concentric shell. A known 
"^quantity of heat per unit time generated electrically 
in the ceramic core flowed radially through the speci- 
men to the concentric ceramic shell. The ceramic 



core had a concentric ring of equally spaced jilioles 
at a radius, r', parallel to the axis, each containing 
a heater wire. Temperatures were measured by an 
axial thermocouple in the ceramic core and by 
tangential tliermocouples in the ceramic shell at 
radius c, onlv one of which is shown in fio-ure 1. 
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Figure 1. Horizontal cross section of the essential elements of 
the apparatus. 

In figure 1, if the cylindrical surfaces r=a and r=b 
are each isothermal, the heat flow between these 
surfaces will be radial except near the ends. In gen- 
eral, the thermal conductivity of the specimen ma- 
terial will vary with temperature, and the heat flow 
rate, q^ through a cylindrical element of unit axial 
length is 



q= — 27rkr 



dr 



(1) 
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Temperature is denoted by the symbol T, the tem- 
perature-dependent thermal conductivity by k, and 
the radial coordinate by r. For reasonable tem- 
perature diflFerences, the thermal conductivity of the 
specimen material can be assumed to vary linearly 
with temperature; then (1) becomes 



- — dT 



(2) 



k is the thermal coiiductivity at an arbitrary refer- 
ence temperature, T, and a is the corresponding 
temperature coefficient of thermal conductivity. 
Integration of (2) yields, after substitution of 
boundary conditions, 



2Tk (T-T,) ^ 
In hja 



(3) 



where Ta is the temperature at r=a, T^ is the tem- 
perature at r=6, and k corresponds to the temperature 
T^ which is given by 

T^_L±^, (4) 

If the circle of heater wires at r=r^ were a contin- 
uous cylindrical heat source, the entire region r<^y 
inside the heater circle would be isothermal. For 
the case of a finite number of line heat sources at 
r=^r\ the temperature measured at the axis is equal 
to the average temperature at the radius of the 
heater circle [1]^ (see also sec. 5). The heat flow 
to the surface of the core per unit length is given by 



^2arkJJ\-T^ 



(5) 



where To is the temperature measured in the axial 
well, and kc is the thermal conductivity of the 
ceramic core material, corresponding to a tempera- 
ture (ro+TJ/2. 

Similarly, for the outer cenunic shell, 



2irk,(T,-Tc) 
' lnc/6 ' 



(6) 



where Tc is the temperature measured at r=c, and 
kg is the thermal conductivity of the ceramic shell 
material at (Tft+T,)/2. 

If there are m heater wires at r=r\ each generating 
heat at the rate (/ per unit length, the total heat 
generated per unit length of core is 

g=mg^ (7) 

Combination of eqs (3), (5), (6), and (7) yields 

Y__ "fnq^ In b/a 



2. |r„-r.-— (^^-^+^1 



(8) 



1 Figures in brackets indicate the literature references at the end of tliis paper. 



which is the equation used for purposes of calcula- 
tion of all thermal conductivity data presented in 
this paper. 

3. Experimental Procedure 
3.1. Apparatus 

A vertical cross section of the thermal conductivity 
apparatus is shown in figure 2. The specimen was 
contained in the annular space between the central 
ceramic core, which was supported at both ends, and 
the concentric ceramic shell, which was supported by 
Transite end pieces. The space between the ceramic 
shell and the outer stainless steel case contained a 
loose-fill thermal insulation. 

The ceramic shell was a hollow mulHte cylinder, 
61 cm long, with an inner radius of 3.2 cm and an 
outer radius of 3.8 cm. A helical groove (about 1.6 
turns/cm) on the outer surface was wound with three \ 
separate heaters of 0.10 cm Kanthal A-1 heater 
wire. The main shell heater was 41 cm in length; 
the upper and lower shell heaters were 5 cm in 
length. A length of 5 cm at each end was not 
wound. The upper and lower shell heaters pro- 
vided extra heat needed to offset end losses, so that 
longitudinal temperature differences in the central 
region could be minimized. 

The central core was an extruded mulhte rod, 46 
cm long and 1.27 cm in diameter. Sixteen equally 
spaced lioles, 0.08 cm in diameter, extended the 
entire length of the rod. The centers of the holes 
formed a circle of 0.43 cm radius. The core heater, 
which provided the heat flowing radially through the 
specimen, consisted of a continuous length of plati- 
num-20 percent rhodium (0.05 cm diam) wire 
threaded back and forth through the 16 holes. 
Current lead wires to this heater, and lead wires to 
its two voltage taps located at the upper end of the 
ceramic core, were brought in through the hollow 
ceramic upper support for the core. The core was 
supported below by a ceramic pin, as shown in 
figure 2. Both the upper and lower supports were 
provided with small platinuni-20 percent rhodium 
heaters to minimize longitudinal heat flow at the 
ends of the ceramic core. 

Temperatures in the apparatus were determined by 
means of six platinum: platinum-10 percent rhodium 
thermocouples. These thermocouples were fabri- 
cated from calibrated thermocouple wire, 0.04 cm 
in diameter. 

Three thermocouples were positioned in the 
ceramic shell at the midplane of the apparatus. 
These Avere inserted into tangential holes (fig. 1) 
located 120° apart at a radius of 3.5 cm. The 
temperature, T^, in (6) and (8) was obtained by 
taking the average of the temperatures indicated 
by these three thermocouples. Two additional «■ 
tangential thermocouples were located in the ceramic 
shell, one 15 cm above and one 15 cm below the 
midplane of the apparatus. 

Temperatures in the core were measured by means 
of a thermocouple which was located in a 0.25-cm 
diam axial well and which could be moved vertically 
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Figure 2. Vertical cross section of the apparatus for deter- 
mining the thermal conductivity of powders. 



A. Terminal head for current 

voltage leads. 

B. Ceramic support tube. 

C. Upper support heater. 

D. Upper shell heater. 

E. Main shell beater. 

F. Ceramic core. 

G. Specimen, 



and II. Ceramic shell. 

I. Shell insulation. 
J. Stainless steel case. 
K. Lower support heater. 
L. Lower shell heater. 
M. Ceramic support rod. 
N. Removal)le plug. 
T. Thermocouple. 



by exterior iiianipubition. A scale fixed above the 
apparatus indicated the position of the thermo- 
couple junction. The temperature measured at 
midlength was desii2:nated as Tq. 

3.2. Instrumentation 

The three heaters on the ceramic shell were fed by 
variable voltage transformers, which in turn were ^^ 
fed by a voltage regulating transformer. These 
heaters were individually regulated by on-off con- 
trollers actuated by thermocouples located in the 
ceramic shell adjacent to the heater windings. The 
heaters on the ceramic core supports were manually 
adjusted by means of variable voltage transformers. 

The ceramic core heater, which provided the heat 
flowing radially through the specimen, was powered 
by a 28'V, 4-amp, regulated d-c power supply, 
with a small bank of power resistors in series for 
adjusting to the desired heater current. Power 
input to the ceramic core heater was determined by 
measuring the d-c current through the heater and 
the voltage (h'op across the entire length of lieater 
wire in tlie core. These measurements were made 
using calibrated shunt and volt boxes and measur- 
ing their output voltages by means of a precision 
d-c potentiometer. 

The noble metal leads of the six measuring thermo- 
couples were brought to an isothermal zone box at 
room tempei-ature. A thermocouple with one junc- 
tion in tlie zone box and one in an ice bath was 
placed in series with a, double-pole selector switch, 
so that each measuring thermocouple was auto- 
matically referenced against the ice bath |21. Ther- 
mocouple voltages were read to 0.1 fxv on the preci- 
sion potentiometer. 

3.3. Test Procedure 

The specimen material was placed in the annular 
sj)ace between the ceramic core and the ceramic shell 
til rough the opening at the top of the apparatus. 
Tlie bulk density of the specimen in place was com- 
puted fj'om the weight of specimen material and the 
Vvolume of the test chamber. 

The shell temperature controllers were set at the 
desired temperature, which for the first test in any 
series was about TOO °C, and all heaters were turned 
on. The controllers and power for the heaters at 
the top and bottom of the ceramic shell were ad- 
justed, if necessary, until the temperatures indicated 
by the upper and lower tangential thermocouples 
agreed to within 1 or 2 deg C with the average 
temperature of the three midplane tangential thermo- 
couples. The power to the core heater was adjusted 
to attain the desired radial temperature drop across 
tlie specimen, and the power to the heaters on the 
ceramic core supports was adjusted until a traverse 
of the thermocouple in the axial hole indicated the 
desired longitudinal temperature distribution. 

When a satisfactory steady-state condition was 
obtained, the test data were recorded. With the 
axial thermocouple in the midplane position, three 
sets of readings (at about 20 min intervals) were 
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taken of all tlierniocouples and of the output of the 
volt box and shunt box. At the end of these read- 
ings, the axial thermocouple was moved to a position 
near the top of the ceramic core. After this thermo- 
couple again reached thermal equilibrium, its output 
was read on the potentiometer. Readings were 
taken at 5 cm intervals over a region 15 cm or more 
on either side of the midplane. When a complete 
set of data had been taken, the apparatus was 
adjusted to a new temperature level, as desired for 
the next test. 

3.4. Calculation of Results 

Calculations of thermal conductivity values by 
means of eq (8) were performed by means of a 
digital computer. The data input to the computer 
consisted of the average emf readings for each 
thermocouple and the average readings from the volt 
and shunt boxes. Equations for the estimated 
thermal conductivities of the ceramic core and shell 
as functions of temperature, and for temperatures 
as a function of thermocouple emfs, were contained 
in the computer progi'am. 



4. Results 
4.1. Diatomaceous Earth 

Test results for a specimen of commercial dia- 
tomaceous earth are presented in figure 3. The 
in-place bulk density of the specimen was 0.15 g/cm^. 
The temperature difference across the specimen 
during the various tests was from 30 to 60 deg C 

The open squares in figure 3 represent results 
obtained in five tests before heaters were placed at 
the ends of the ceramic core to minimize axial heat 
flows. For these five tests, the temperature differ- 
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Figure 3. Thermal conductivity of a diatomaceous earth as a 
function of temperature. 



ence between the ceramic core and shell was 10 to 
25 percent less at positions 15 cm from the mid- 
plane than it was at the midplane. (These data, as 
well as those presented below, were calculated using 
eq (8) ; the effect of a nonuniform axial temperature 
distribution is discussed in sec. 5.2.) 

After these five tests, the specimen was removed 
and heaters were installed on the ceramic core 
supports. The same sample was then reinstalled 
as the specunen, and the data represented by the 
open circles were taken in order of increasing mean 
temperature. After the test at a mean temperature 
of 946 °(\ tests were conducted at several lower 
temperatures and finally at a maximum temperature. 
The results of these tests, indicated in figure 3 by 
circles with vertical lines through them, were sig- 
nificantly higher than those of previous tests. 
Removal of the specunen indicated that a slight 
sintering had occurred.^ The powder appeared to 
be slightly fused together, but was easily broken 
apart. The same sample was pulverized by hand, 
replaced in the apparatus, and the test results indi- 
cated by the open triangles were obtained. 

The curve drawn through the data points in 
figure 3 is the quadratic equation of least-mean- 
squares fit to the open squares, circles, and triangles. 
The equation for this line is 



A:=0.414 + 1.114 



(w)+o-o^KmT 



(9) 



where k is expressed in milliwatts/cin-deg C and T 
in °C. For all data points except the circles with 
vertical lines through them, the standard deviation 
of a point, estimated from the residuals about the 
fitted line, is 0.013 mw/cm-C, corresponding to 
1.3 percent at 500 °C. 

The solid ciicles shown in the lower left hand 
corner of figure 3 represent results of tests made on 
ii sample of diatomaceous earth from the same 
container, at the same density, in the NBS guarded 
hot-plate apparatus (ASTM C177). 

Figure 4 presents some literature data on dia- 
tomaceous earth of various densities [3, 4]; the curve 
from figure 3 is reproduced for comparison. There 
is substantial agreement between tlie results of 
this investigation and those from the literature, 
taking density into consideration. 

4.2. Powdered Alumina 

Test results for an aluminum oxide powder are 
presented in figure 5. This alumina w^as produced 
at the National Bureau of Standards by ignition 
of hydrated aluminum chloride in a muffle furnace 
at 1,150 °C. The method of production and some 
of the properties of this type of alumina powder 
have been described [5]. 



2 An uncontaminated diatomaceous earth should melt at the temperature of 
fusion of silica, which is well above any temperatures obtained in these tests. 
A commercial diatomaceous earth, however, may contain sufficient clay to 
produce sintering at temperatures as low as 800 °C. No tests were conducted 
to determine the actual temperature at which sintering began for the diatomaceous 
earth discussed in this paper. In general, diatomaceous earth should be used 
with caution at temperatures above^SOO °C if sintering would cause difficulties. 
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INT. CRITICAL TABLES 
CALVERT AND CALDWELL 
NBS SMOOTHED DATA 
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Figure 4. 



Comparison of diatomaceoiis earth thermal conduc- 
tivity results with literature data. 
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Thcrmcd conductivity of an aluminum oxide poivder 
as a function of temperature. 



The data indicated by the open squares in figure 5 
were taken in order of increasing mean temperature 
on a specimen having a bulk density^ of 0.44 g/cm^. 
The curve drawn through the squares is the quadratic 
equation of least-mean-squares fit. The equation 
for this line is 



= 1.130+0.566 (^)+0.333 



T Y 

1000/ ' 



(10) 



where k is expressed in mw/cm-C and T in °C. The 
estimated standard deviation of a point is 0.010 
mw/cm-C, corresponding to 0.7 percent at 500 °C. 



3 Later clucks (if tlio packing procoduro used for this specimen indicate a pos- 
sible uiiccriiii mi. \ (if 0.02 u/cHi'* ill the rei)()rte(i iii-place bulk density. All other 
densities irpornd arc believed to be aecunite lo within 0.005 g/cm3. 



The data indicated by the open circles were taken 
in order of increasing mean temperature on a speci- 
men having a bulk density of 0.40 g/cnr'^. The speci- 
men was the same sample as that previously tested, 
but w^as intentionally packed to a lower bulk density. 
After tlie tests up to 919 °C, a test, indicated by an 
open triangle, was made at 526 °C. The curve drawn 
through the open circles and the triangle is the 
quadratic equation of least-mean-squares fit to these 
points. This equation is 



k=0M7 +0.841 (i^)+0-2% ( 



1000 



)' 



(11) 



in the same units as above. The estimated standard 
deviation of a point is 0.011 mw/cm-C, corresponding 
to 0.8 percent at 500 °C. 

The solid circles shown in the lower left-hand corner 
of figure 5 represent results of two tests made in the 
NBS guarded liot-plate apparatus on a sample of 
alumina pow(k>r, taken from the same batch as the 
sample installed in the ceramic core apparatus, liav- 
ing a bulk density of 0.41 g/cnr^ The solid diamond 
shown represents a guarded hot-plate test on the 
sample of alumina which liad been tested in the 
ceramic core a})paratus at higher temperatures. The 
in-|)lace l)ulk (hnisity for this test was 0.46 g/cm^. 

The higher density specimen (0.44 g/cm^) had a 
room temperature thermal conductivity about 20 
percent liigher than that of the lower density speci- 
men (0.40 g/cm^). However, the curves corresj)ond- 
ing to these different densities converged with in- 
creasing temperature and crossed at about 730 °C. 
This crossing of the curves is believed to arise from 
increased radiation within the larger pores of the 
lower density material. 

5. Discussion of Test Method 

5.1. Determination of Temperature Difference Across 
the Specimen 

Measurement of the temperature diflference across 
the specimen, using the core and shell thermocouples 
shown in figure 1, represents a departure from the 
procedure used in radial flow methods for loose-fill 
materials described in the recent literature [6, 7]. 
The possibility of chemical contamination of thermo- 
couples by the specimen material is greatly reduced, 
and the often-difficult problem of determining and 
maintaining the distance between thermocouples lo- 
cated within the specimen material is avoided. How- 
ever, with the method of measurement used here, 
account must be taken of temperature drops in the 
core and shell, and at their interfaces with the speci- 
men if necessary, to obtain tlie temperature drop 
across the specimen. 

The thermocouple in the axial hole at midlength 
of the central ceramic core is located in a region 
where there are substantially no temperature gra- 
dients, and therefore the thermocouple accurately 
attains the temperature of its environment without 
disturbing the heat flow pattern. The theory of a 
harmonic logarithmic potential function shows that 
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the value of the temperature at the center of a circle 
is equal to the average of its values over the circum- 
ference of the circle [8,9]. Hence, the temperature 
measured at the axis of the core is equal to the 
average temperature at the fictitious surface, r=^r' . 
The difference between the average temperature of 
the fictitious surface at r=r^ and the average tem- 
perature of the surface of the core at r—a may be 
computed from (5). It is to be noted that eq (5) is 
valid regardless of the distribution of heat sources 
on the circle T=r\ as has also been shown by Peavy 
[1] in analyzing the two-body system of the specimen 
and the core with line heat sources. 

Similarly, the temperature difference, Tf^—Tc, 
from the inner surface of the shell to the radius of 
the tan2:ential thermocouples may be computed 
from (6)^ 

The uncertainty in determining the temperature 
difference between the axial thermocouple in the 
core and the three thermocouples in the shell is 
estimated"^ to be less than 1 percent. The computed 
temperature drops in the core (Tq—To) and in the 
shell (Tjj—Tc) were less than 2.4 percent and 0.6 
percent, respectively, of the temperature drop 
through the specimen for all of the results presented 
in section 4, using the best available values for the 
thermal conductivities of the mullite core and shell, 
the uncertainties in which may have been as much as 
30 percent. 

The effective radius, r\ of the virtual heating 
surface in the core is uncertain by not more than the 
radius of the heater holes, thus introducing an addi- 
tional error of not more than 25 percent of Tq—Tq. 
Similarly, the effective radius, c, to the tangential 
thermocouple wells is uncertain by not more than 
the radius of the wells, thus introducing an additional 
error of not more than 35 percent of T^—T^. Com- 
bining the errors arising from the uncertainties in kc 
and ks and those from In a/r^ and In c/6, the esti- 
mated uncertainties in T^—Ta and T^—Tc are 40 
and 50 percent, respectiveh', thus introducing an 
additional uncertainty in Ta~ T^ of 1.3 percent. For 
fine powders of low thermal conductivity, such as 
those tested, errors due to thermal contact resistance 
between the specimen and the surfaces of the core 
and shell are believed to be negligible. Thus, com- 
bining uncertainties, the temperature drop across the 
specimen (the expression in braces in the denominator 
of eq (8)) is considered to be known with an uncer- 
tainty of 1.5 percent. 

5.2. Longitudinal Heat Flow 

Equation (8) was derived on the assumption of no 
longitudinal heat fiow\ For an apparatus of finite 
length, having imperfect end guarding, the effects of 
axial heat flow should be considered. Jakob [11] 
considered the effect of axial heat flow for the sym- 
metrical case in which both ends of the core are at 
the same temperature but are at a temperature 
different from that at the center. An analysis is 
presented here which allows different temperatures 

4 All ostimatod uncortainties in this paper represent the author's estimate of 
the "two-signia" limit. All uncertainties are combined using • tandard propaga- 



tion of error formulas (see, for instance, Davies [10]). 



at the two ends of the core, and also takes into 
account the effect of core temperature distribution 
on the power measurement if the electrical resistiv- 
ity of the heater wire varies linearly with temperature. 
A simplified cross-sectional viev*' of the apparatus 
is shown in figure 6. It is assumed that the surface 
r=b is isothermal at temperature T^ and that tem- 
peratures are measured in the core at three positions 
equally spaced along the axis. It is assumed that, 
for the purpose of analyzing longitudinal heat flow, 
a transverse cross section of the core can be treated 
as being isothermal. 



SPECIMEN 

\ 



= 0l 

Z = -M 



Figure 6. General core temperature profile used in analyzing 
effects of longitudinal core temperature imbalance. 



Neglecting longitudinal heat flow in the specimen, 
the differential equation for heat flow along the core 
can be written as 



(Pv 27rkv 



mq 



dz' C In b/a 



C 



(12) 



where 



v=T—J\, T being the temperature at any 

point along the core; 
s=loiio-itudinal coordinate ; 



and 



C=the longitudinal thermal conductance per 
unit length of the core, including heater 
wires and holes. 

In general, the rate of heat generation, g', per unit 
length of heater wire will not be constant if the elec- 
trical resistivity of the heater wire varies appreciably 
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with temperature. For reasonable tempeiature var- 
iations, the electrical resistance ])er unit leng'tli of 
heater wire can be assumed to l)e a linear function 
of temperature, 

P=Po{l+yv), (13) 

where p and po are the electrical resistances of a unit 
length of heater wire at temperatures T and Tf„ re- 
spectively, and 7 is the fractional change in electrical 
resistivity per unit change in temperature. The rate 
of heat generation per unit length of heater wire is 



(/=Pp=Ppo(l+yv), 



(14) 



where I is the electrical current flowing through the 
heater. Equation (14) must be substituted into (1 2) . 
The conditions which are to be satisfied are 



v=Vi at z^ 
v=Vq at z--^ 
v=V2 at z^ 



Substitution of these conditions 
(12) yields 



(15) 
into tlie solution of 



^=Fo+(- 



ill wliicli 



F,+F2-2Fo\e,osliM.--l 



) 



cosh jiL 



■1 



\ 2 / Sinn 11 L 



^-^[ 



2irk 



The thermal 
given by 

7_ mppojn^i/o) 



In {hi a) 
conductivitv 



ym]' 



Oi) . 



(17) 



of the s[)(M*ini(Mi is 



[ 



1+yFo- 



T',+F,-2F„ 



2FoC.s1.m/.-(Fi + 



(18) 



The total power delivered by the core heater is 
determined by the current, /, through the heater 
and the voltage drop, E, across the total length of 
heater wire in the core. The total resistance of the 
heater is given b}^ 



R=E!I=p,{l+yV)S 



(19) 



where S is the total length of heater wire and V is 
the average temperature of the core heater, obtained 
by integration of (16) over the entire length, 2M, 
of the core, 



— 1 C^^ 



+ 



■2Fo \ sinh mM 
/m-M(cos1i nL—l) 



(20) 



k^k 



Eciuiition (18) can be rewritten as 

F,+F2-2F„ r 1 



. fj , F1+F2-2F0 r 
L 1+7F I 



2Fo cosh Mi- (F.+Fi) 
7(sinli fxM—ixM) 



2MM(cosh/iL— 1) 



}] 



(21) 



where 



j,_ mEI\n (b/a) 
" 27rFo'S 



(22) 



is the apparent value of thermal conductivity ob- 
tained if no corrections are made for the effects of 
axial heat flow. Since k is defined in terms of /x and /x 
in terms of k, iterations are necessary in order to 
obtain k. A first approximation to k is given by k^; 
substitution of k^ into (17) yields a value of ijl wliich 
can be used in (21). The resultant k can be used to 
obtain a better value of ju, etc. The convergence is 
quite rapid in all practical cases; the first approxi- 
mation did not introduce an inaccuracy greater than 
0.001] ])ercent in k foi* the case to be discussed below. 
E(iuation (22) corresponds to (8), the equation that 
was used for purposes of calculation; EI/S corre- 
sponds to (j\ and V'o corresponds to the expression 
in braces in the denominator of (S). The first expres- 
sion in the braces in (21) represents the correction 
for axial iieat flow; the second expression represents 
tlie correction in power measurement to take account 
of the variation of electrical resistarux> with tem- 
perature. It is of interest to note that, for a ])Ositivc 
value of 7, these corrections are of oi)posite sign. 

Data were presented in section 4.1 for which the 
temperature difl'eren<'e between tlie ceramic core and 
the substantially isothermal ceramic shell was 10 to 
25 percent less at positions 15 cm from the midplane 
than at the center. As seen in figure 3, the results 
for these tests were not appreciably differeiit from 
the conductivities obtained later when the core w^as 
substantially uniform in temperature along its length. 

In order to investigate more thoroughly tlie effects 
of axial heat flow^ in the apparatus, a series of five 
tests was conducted in w^hich the ends of the ceramic 
core were held at temperatures different from that in 
the core at the midplane. These tests were all con- 
ducted at a mean temperature of 525 °C, with the 
alumina powder specimen discussed in section 4.2. 

The ceramic shell was maintained approximately 
isothermal over a central region of at least 30 cm 
length. The core temperature at positions 15 cm 
above and below the midplane was varied by adjust- 
ment of the smrdl heaters at each end of the ceramic 
core. The results of these tests are shown in figure 
7, where the percent departure of the thermal con- 
ductivities from their mean value is plotted against 
the quantity (Fi+ F2)/2yo. The solid circles repre- 
sent experimental data; the solid curve is the straight 
line of least-mean-squares fit to the data. A straight 
line w^as used because, for the various values of pa- 
rameters involved, tlie departure can be shown to be 
a nearly linear function of iVi+V2)l2Vo for values 
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not too different from 1.0. The slope of tlie straight 
line in figure 7 was used to compute the value for ju 
of 0.29 cm"^ where 7 was taken as 8.0 X 10"^ deg 
C\ Vq as 44 deg C, L as 15 cm, and M as 23 cm. 
Substitution of this value of m into (17) yields for 
the longitudinal conductance of the core, C, a value 
of 0.06 w-cm/C. This value agrees with the esti- 
mated longitudinal conductance of tlie ceramic core 
and heater wires. 






(V,+ V2)/2Vo 

Figure 7. Effect of longitudinal core temperature imbalance, 
(Vi + V2)/2Vo, on thermal conductivity values obtained for 
alumina powder. 

With the exception of those tests in which core end 
heaters were not used, all data presented in section 
4 gave values of (Fi+ V2)/''^Vo between 0.85 and 1.15. 
From figure 7, it is seen that this corresponds to an 
uncertaint}^ of less than 0.4 percent in the thermal 
conductivity values obtained using (8). Values of 
(Vi+V2)/2Vc, could be held between 0.95 and 1.05 
without undue difficulty, if desired, thus further 
reducing this uncertainty. 

5.3. Eccentricity 

It has been assumed that the outer surface of the 
ceramic core and the inner surface of the ceramic shell 
were concentric cylinders. Since concentricity may 
not always be possible, the effect of eccentricities 
should be considered. The thermal resistance be- 
tween two eccentric parallel circular cylinders is 
given by Carslaw and Jaeger [12] and is, in the 
notation of this paper, 



T.-T, 



= — - arccosh 



2Trk 



2ab 



(23) 



where e is the eccentricity, or the distance between 
the centers of the two cylinders. 

From (23), the thermal conductivity is given by 



k= 



MTa-T,) 



arccosh 



l+X^-e^(X-l)^ 
2X 



(24) 



where, following El-Saden [13], the dimensionless 
eccentricity e=e/{b—a) varies between zero (con- 
centricity) and one (cylinders in contact), and 
\=r-b!a is the ratio of the radii.^^ 

For the case of eccentric cylinders, if (3) is used 
instead of (24) to compute the thermal conductivity 
of the specimen, an error will be introduced. Since 
it is of general interest, the percent error due to 
neglecting the effect of eccentricity is shown in 
figure 8 plotted against dimensionless eccentricity, e^ 
for several radius ratios, \^b/a. 




Figure 8. (Calculated error in thermal conductivity due to 
neglecting the effect of eccentricity for radial heat flow between 
isothermal cylinders. 

For the apparatus described in this paper, the 
ratio X was approximately 5.0. The dimensionless 
eccentricity as defined above was less than 0.05, so 
that the uncertainty due to eccentricity was less 
than 0.1 percent. 

5.5. Precision and A'^tTiracy of Resultti 

The precision of measurement is indicated by the 
standard deviations (estimated) of the data, which 
are of the order of 1 percent. 

The accuracy of the reported measurements is 
more difficult to estimate. The total power gener- 
ated in the core heater was known to within 0.1 
percent. The length of wire between the potential 
taps at room temperature was known to within 0.2 
percent. Thermal expansion could have increased 
this length b}^ 1 percent at 1,000 °C if friction 
between the heater wire and the ceramic core at 
points of contact did not prevent free longitudinal 
expansion of the wire in the holes. Because the 



fi Equation (24) can also be derived from the last equation presented by El- 
Saden. To convert El-Saden's eq (38) to the form presented abo ve, it is helpf ul 
to make use of the identity arccosh r -arccosh y= arcco sh (xy— V(a:2— 1)(2/2— 1)). 

« Making use of the identity arccosh a: = In Ct+V^:^— 1), eq (24) reduces to the 
equivalent of (3) for the case of concentric cylinders (€=0). 
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effect of oxpaiisiou was uncertain, no expansion 
corrections were made. The heat flowing- radially 
throug'h the specimen is i'm'tlier in question by not 
more than 0.4 percent in an extreme case, due to the 
uncorrected effects of longitudinal temperature vari- 
ations along the core (see sec. 5.2). 

In section 5.1, it was estimated that temperature 
drops across the specimen were determined with an 
uncertainty of 1.5 percent. Uncertainties in the 
ratio of the sliell radius to that of the core did not 
introduce more than 0.2 percent uncertainty in the 
results. The error due to possible eccentricity of 
the ceramic core in the ceramic shell did not exceed 
0.1 percent. Combining all of the uncertainties 
discussed, uncertainties in the thermal conductivity 
results presented in section 4 are estimated to be 
less than 2 percent at 100 °C, and, because of heater 
wire expansion, 3 percent at 1,000 °C. 

5.5. Comments 

It is felt that the type of central heater used in 
this apparatus has advantages which juay recom- 
mend it to others interested in radial heat flow 
measurements. The particular advantages of this 
arrangement are, (1) the measuring thermocouple in 
the core attains the temperature to be measured and 
does not interfere with the pattern of heat flow, (2) 
the measuring thermocouple is well protected both 
mechanically and clieniically but still can be easily 
removed for calibration or renewal, and (3) tem- 
peratures can be measured at any desired longitudinal 
position using the same thermocouple. 

By means of suitable modifications, the e<[uip- 
ment described could be extended to cryogenic 
temperatures, or to more elevated temperatures. 
If required, the accuracy of measurement could be 
improved over that attained in the apparatus de- 
scribed by taking careful account of all the factors 
discussed in sections 5.1 and 5.2. Api)lications for 
which a heated core of this type might be advan- 
tageous include determination of forced or natural 
convection from cylinders and deternunation of the 
thermal conductivity of pipe insulations or fluids. 
In the latter case, convective effects would need 
to be excluded by extra|jolating results obtained with 
small temperature differences to zero temperature 
difference. 



Tlie experimental measurements were performed 
by C. I. Siu and J. E. Griffith. The autlior thanks 
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